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During the course of a gravitational lensing survey of distant, X-ray selected, EMSS clusters 
of galaxies, we have studied 6 X-ray luminous (L x > 5xl0 44 h^Q erg s _1 ) clusters at redshifts 
exceeding z = 0.5. All of these clusters are apparently massive. In addition to their high X- 
ray luminosity, two of the clusters at ,?~0.6 exhibit gravitationally lensed arcs. Furthermore, 
<3\ the highest redshift cluster in our sample, MS 1054—0321 at z = 0.826, is both extremely X-ray 

luminous (£o.3-3.5keV = 9.3x 10 44 h^Q erg s _1 ) and exceedingly rich with an optical richness 
comparable to an Abed Richness Class 4 cluster. 

In this Letter, we discuss the cosmological implications of the very existence of these clusters 
for hierarchical stucture formation theories such as standard = 1 CDM, hybrid = 1 C+HDM, 

r-5 and flat, low-density A + CDM models. 
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Subject headings: cosmology: observations — large-scale structure of the universe — 

^ gravitational lensing — galaxies: clusters: individual (MS 1054—0321) — 

Ci X-rays: galaxies 



1. Introduction 

The very existence of massive clusters at high redshift (z > 0.5) is problematic for standard cold dark 
matter (CDM) theories of hierarchical structure formation. The CDM model, as originally formulated, 
assumed that the large scale distribution of galaxies could be used to trace the underlying mass. It became 
clear early on, however, that observed structure and galaxy velocities on small scales demanded a biased 
model of galaxy formation where the galaxy number fluctuations SN/N are more pronounced than the 
actual mass fluctuations SM/M by some biasing factor b > 1 (Davis & Peebles 1983). This bias factor is 
defined as the ratio of the galaxy to mass fluctuations within an 8 h~ x Mpc radius sphere. Since the field 
galaxy number fluctuations on this scale are observed to be close to one (Davis & Peebles 1983), the bias 
parameter is often expressed as b=l/ag, the inverse of the rms density fluctuation spectrum filtered on 
this 8 h~ x Mpc scale. Large clusters of galaxies are rare objects that formed from exceedingly large mass 
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density perturbations. Consequently, the abundance of clusters as a function of mass provides a sensitive, 
and nearly model free way to measure the CDM spectrum normalization. A factor of two change in <jg 
alters the expected number density of clusters by two orders of magnitude (Carlberg et al. 1993; Bahcall 
& Cen 1992; White et al. 1993). Furthermore, the change in the comoving number density of rich clusters 
as a function of redshift provides another important observational constraint on the normalization of the 
CDM spectrum. The number density of clusters evolves faster in a high bias CDM universe than in a low 
bias CDM universe (Cen k Ostriker 1994). 

Various groups have attempted to measure and analyze the abundance of rich clusters of galaxies at both 
high and low redshift. Bahcall & Cen (1992) find that cluster abundances at low- z are inconsistent with 
= 1 CDM for any value of <jg, and are best fitted by <jg ~ 1, low density (Oq ~ 0.25) CDM models, either 
with or without a cosmological constant. Conversely, White et al. (1993) argue that observed low- z cluster 
abundances suggest <jg ~ 0.6 for an Oq = 1 universe. At higher redshift, cluster abundances have been more 
difficult to determine, but results seem to indicate an overabundance of massive clusters, again in favor of 
a low, or even unbiased CDM model with <jg ~ b ~ 1. For example, Evrard (1989) argued that the existence 
of only three clusters with high velocity dispersions in the intermediate redshift range of z ~ 0.4 — 0.6 is 
difficult to reconcile with a high value for b. Peebles, Daly & Juszkiewicz (1989) showed that the existence 
of just one cluster of ~ 10 15 Mq in the redshift range z ~ 0.5 — 0.6 (e.g. CL 0016+16) implied an upper limit 
of b < 1.6. Their limit drops to b < 1.38 if just one Coma-like cluster is found between redshifts z ~ 0.7 — 0.8. 
And such clusters do appear to exist. Gunn (1990) lists 4 clusters in the range 0.7 < z < 0.85 and one at 
z ~ 0.9, two of which have measured velocity dispersions of ~ 1000 km s _1 (i.e. comparable to Coma). 

Of course the statistics on the abundance of high- z clusters are poor because the presently-known num- 
ber of such clusters is so small. Furthermore, one might argue that many of these clusters may not be as 
massive as they appear due to inflation of the measured velocity dispersion from field galaxy contaminaton 
and projection effects (Frenk et al. 1990). Recent optical and X-ray observations of some optically-selected 
clusters with z > 0.6 may lend support to these contamination arguments. Thimm & Belloni (1994) pre- 
sented observations of two of the highest redshift Gunn, Hoessel & Oke (1986) clusters and found that one 
of them (CL 2155+0334) was not a cluster at all, but a chance superposition of field galaxies. Castander et 
al. (1994) observed 10 distant, optically-selected clusters (from Gunn 1990) with ROSAT. Only 5 of their 
clusters had known redshifts, but all of those 5 were at z>0.7. Castander et al. detected X-ray emission 
from only 2 (and possibly a third) of these clusters. After converting their detected fluxes to luminosities in 
the standard way using aT~7 keV thermal spectrum and correcting for the extended emission assuming 
a (3 = 2/3 model for the cluster density profile, Castander et al. found the detected clusters were relatively 
modest X-ray emitters with much lower X-ray luminosities (£o.i-2.4keV ~ 1 X 10 44 h^Q erg s _1 ) than 
equivalently-rich low or intermediate redshift clusters. In a similar study, Nichol et al. (1994) conducted a 
ROSAT search for distant X-ray clusters in fields containing either previously-identified, optically-selected 
clusters or wide angle radio sources. Two of their objects were detected with ROSAT; one corresponded 
to one of their optically-selected clusters (z=0.61), while the other corresponded to a radio-selected cluster 
estimated (from photometry) to be at z~0.7. Again, after converting the detected fluxes to luminosities 
in approximately the same manner as Castander et al. , Nichol et al. found that both of their clusters 
were also relatively modest X-ray emitters (T x ~ 10 44 erg s _1 ). Nichol et al. estimated a lower limit to 
the comoving number density of such clusters (n z > 1.2xl0 -7 h\ Q Mpc -3 ) and argued this was marginally 
consistent with a standard b ~ 2 CDM universe. 

While these data would seem to imply that extremely rich systems with large amounts of X-ray emitting 
gas do not exist at these high redshifts, we will show several examples that indicate this is not the case. 
In this paper, we present preliminary data on 6 rich, massive, X-ray luminous clusters at z > 0.5 (z = 0.66, 
z max = 0.826). The clusters were selected from the Einstein Observatory Extended Medium Sensitivity 
Survey (EMSS; Gioia et al. 1990a) and are the most distant members of a sample of 40 such clusters we are 
examining for evidence of gravitational lensing (Luppino & Gioia 1992; Luppino et al. 1995). In the following 
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section, we briefly discuss the EMSS cluster sample, and describe some of the properties of the six z > 0.5 
clusters. In §3, we discuss the implications the very existence of these clusters have for various CDM models. 
Unless specifically stated otherwise, all relevent cosmological quantities were computed assuming Ho = 50 
and = 1, however we will often list quantities in terms of h^o = Ho/ 50 to show the explicit dependence of 
that quantity on the value of the Hubble parameter. 

2. Highest redshift EMSS clusters 

The EMSS cluster sample consists of 103 clusters of galaxies extracted from a total of 835 EMSS sources. 
The most recent optical, X-ray and radio data for these clusters have been compiled and presented by Gioia 
& Luppino (1994). The cluster subsample chosen for our gravitational lensing survey (Luppino et al. 1995) 
was subjected to the following selection criteria: 1) the fluxes of the sources had to exceed 1.33xl0 -13 erg 
cm -2 s _1 , 2) the sources had to lie North of declination 8 = —40°, and 3) the clusters had to be distant 
(z > 0.15), and 4) X-ray Luminous (L x > 2.0xl0 44 h^ 2 erg s _1 ). Note, X-ray luminosities were determined 
in the conventional way assuming aT~6 keV thermal spectrum and correcting for the extended emission 
assuming a (3 = 2/3 model for the cluster surface brightness (the reader should refer to Gioia & Luppino 
1994 and references therein for the details of this procedure). These selection criteria resulted in 40 clusters 
of galaxies spanning a redshift range of 0.15 < z < 0.83, and included all of the EMSS clusters with z > 0.5. 

Optical imaging observations of these clusters were carried out using the UH 2.2 m, with some data 
obtained using the CFHT and Keck telescopes. The details of the imaging observations and data reduction 
can be found in Luppino et al. (1995). The majority of our spectra were obtained using the MOS on the 
CFHT or LRIS on Keck. Details of the spectroscopic observations and data reduction will be presented in 
a later paper. 



Table 1. z > 0.5 EMSS clusters 
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Koo 1981; b Dressier & Gunn 1992 



In Table 1, we list the six z > 0.5 EMSS clusters. All 6 of these clusters are extremely X-ray luminous 
(-^ ^ 5x 10 44 h^ 2 erg s _1 ). Two, in fact, are among the most X-ray luminous clusters known — at any 
redshift. In each of the clusters in this paper, we are reasonably confident that the X-ray emission comes 
from the cluster and not from an AGN in the field. Only 2 of the 6 cluster fields have a radio source detected 
at 6 cm with the VLA, and in both cases that source is coincident with the brightest cluster galaxy (see 
Gioia & Luppino 1994). Thus any AGN in these fields would have to be radio quiet. Furthermore, during 
the course of the EMSS optical identifications (Stocke et al. 1991), any blue, stellar objects in the fields were 
observed spectroscopically. While there are several cases in the EMSS where AGN are present in clusters, 
and the identification is not certain, this is not the case for the six clusters discussed in this paper. 

In addition to high X-ray luminosity, there is other compelling evidence that these clusters are genuinely 
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massive. Two of the clusters at 2 ~ 0.6 contain large lensed arcs that allow us to make crude, but independent 
estimates of the projected mass in the clusters cores (Luppino et al. 1995). We find central masses of 
~5x 10 14 M Q enclosed within a 237/i^ 1 kpc critical line for MS 0451-0305 (2=0.55) and ~ 1.3 X 10 14 M© 
enclosed within 121/i"^ kpc for MS 2053-0449 (2=0.583). In both cases, we assume the source redshift 
is twice that of the cluster lens, but these masses drop by only 30% if the source is moved to z s ~ 2. We 
can also estimate the cluster central velocity dispersion from the lensing geometry if we make the simple 
assumption that the cluster potentials are isothermal spheres. In this case, the central velocity dispersion 
(<r) is related to the lensing factor (D s / D\ D\ s where D s , Di, and Di s are the angular diameter distances 
to the source, the lens and from the lens to the source, respectively), and distance of the arc to the center 
of the lens (r c ) by a 2 = (c 2 /47r) (D s / Di Di s ) r c . Again assuming the source redshifts are twice the lens 
redshifts, we find central velocity dispersions of ~ 1750 km s _1 and ~ 1200 km s _1 for MS 0451 — 0305 and 
MS 2053 — 0449 respectively. Note that the masses estimated above are only for the central "enclosed" cores. 
The total masses of each of these two clusters must surely be comparable to or greater than that of Coma. 

Moving even further out in redshift, we find the two most distant EMSS clusters, MS 1137+6625 at 
2 = 0.782 and MS 1054-0321 at 2 = 0.826, are also extremely X-ray luminous (L x ~ 10 45 erg s" 1 ), with X- 
ray luminosities nearly an order of magnitude brighter than the optically-selected 2~0.8 clusters observed 
by Castander et al. (1994) or the optical/radio-selected 2~0.6 clusters observed by Nichol et al. (1994). 
Although neither of these 2~0.8 EMSS clusters show any obvious lensed arcs, their optical appearance 
is striking. MS 1137+6625 is rich and compact (see Fig 1 [Plate 1]), and MS 1054—0321 in particular is 
exceedingly rich (see Figs 2a and 2b [Plates 2 and 3]), with a projected spatial distribution of galaxies 
comparable to the richest low and intermediate redshift clusters such as CL 0016+16 (MS 0015.9+1609) or 
CL 0939+4713 (Gunn et al. 1986; Dressier et al. 1994). We measure a central richness (Bahcall 1981) for 
MS 1054—0321 of A0.5 = 50 + 10 (after background subtraction from the periphery of the CCD frame), which 
is significantly richer than Coma (AQ^ ma = 28; Bahcall 1981) and corresponds roughly to an Abell richness 
class 4 cluster. The A0.5 — a correlation (Bahcall 1981) suggests MS 1054—0321 has a velocity dispersion 
of ~ 2000 km s _1 . We have measured a preliminary value for the velocity dispersion (from spectra of 7 
galaxies) of 1643^44 km s _1 . 

3. Discussion 

In the previous section, we demonstrated that several rich, massive and X-ray luminous clusters do exist 
at moderately high redshifts. In order to obtain a quantitative estimate of the comoving number density of 
these clusters, we take the cosmological expression for the comoving volume element (Mathez et al. 1991), 

dV _ ( c \ 3 (q z + (g - l)(yi + 2q z - l)) 2 

and integrate it over the appropriate solid angle and redshift range of interest. Actually performing this 
integral in the case of the EMSS is a complicated process. We use the 1/Vm ax technique as applied by Gioia 
et al. (1990b) and Henry et al. (1992) in their investigation of the EMSS cluster X-ray luminosity function. 
First we consider the maximum redshift, 2 max , at which a cluster could have been seen for a given solid 
angle that has some associated limiting sensitivity. The search volume for a given cluster is then the sum of 
all volumes within a redshift shell that is bounded on the low end by z\ and on the high end by the lesser of 
z 2 or ^max f° r each different sensitivity limit. The total volume for N clusters over the redshift shell z\—*zi 
is then obtained by summing the volumes for each cluster 
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Vtot = 2^1 

i = l 
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22 (Aim ,14 ) (1 + z - vrr^y 



/lim 



where we have set qo = ^ where 22 is a function of the cluster luminosity and the appropriate survey flux 
limit, and where c?Oj-(/ij m ) is the flux limit dependent sky coverage which ranged from 0.01 steradian for 
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Figure 1. CCD image of MS 1137+6625 at z=0.782. This image is an 800x800 subarray from a 3000s R 
exposure and measures 1.4 h7 Mpc x 1.4 h7 Mpc at the cluster redshift. 
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Figure 2a. CCD image of the center of MS 1054-0321 at z=0.826. This image is an 1024x1024 subarray 
from a 14400 s I exposure and measures 1.9 h7 Mpc x 1.9 h7 Mpc at the cluster redshift. 
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Figure 2b. Wide field CCD image of the field containing MS 1054-0321. The image is a full 2048x2048 
CCD frame and covers a field of view of 7'. 5x7'. 5 or 3.7 h^Mpc x 3.7 h^Mpc at the cluster redshift. 
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/ Um = 1.33 X 10" 13 erg cm" 2 s" 1 to 0.22 steradian for / Um = 3.57 X 10" 12 erg cm" 2 s" 1 (see Table 3 
in Henry et al. 1992). We find a volume of 5.7 X 10 8 Mpc 3 for the redshift shell spanned by the 6 
most distant EMSS clusters (^=0.546 — ► 0.826; ~z = 0.66) indicating a comoving number density at 1=0.66 
of n z (L x > 5 x 10 44 erg s _1 ) = 1.1 ± 0.43 X 10~ 8 /i 3 Mpc -3 . Although this number density is lower than 
the density derived by Nichol et al., we point out that the number density we find is for clusters that are 5 
to 10X more X-ray luminous. Moreover, we consider 1 X 10~ 8 /i 3 Mpc -3 to be a conservative lower limit 
for n z (L x > 10 44 ), particularly since, at z > 0.5, the flux-limited EMSS cluster sample will only include the 
most X-ray luminous clusters, and will miss the ~ 10 44 erg s _1 clusters that no doubt exist, but whose flux 
would have fallen below the 1.33 X 10~ 13 erg cm -2 s _1 EMSS detection criterion. If we restrict ourselves 
to just the two highest redshift EMSS clusters at ~z = 0.8, we find a number density n z (L x > 7.5 x 10 44 erg 
s _1 ) = 7 ±5 X 10~ 9 /i 3 Mpc -3 . We remind the reader, with some trepidation, that these number densities 
are based on only a handful of clusters. Unfortunately, this is an unavoidable limitation of our present 
sample, and the large error bars reflect the small number of clusters used. In Figure 3, we plot the two 
densitites we find at 1=0.66 and ^=0.8 along with the three ,?~0.8 clusters of Castander et al. (1994). 
On this plot we also overlay the integrated power law luminosity function for the largest redshift shell from 
the EMSS study of the X-ray luminosity function (Henry et al. 1992). It appears that the X-ray luminosity 
function at these high redshifts is consistent with the LF for EMSS clusters in the 0.30 < z < 0.60 {~z = 0.33) 
shell. 

To compare these number densities to the predicted abundances of clusters for CDM models, we take the 
z ~ mass functions of White et al. (1993) or of Bahcall & Cen (1992), and shift them to high- z for various 
values of <jg following Cen & Ostriker (1994). We argue that all of our z > 0.5 EMSS clusters are Coma-like, 
and thus are of order 10 15 Mq within the central 1 Mpc. As we stated earlier, we believe this argument is 
valid given the estimated core lensing masses in two cases, and the extremely high L x in the two additional, 
highest-,? cases. We find the number density of clusters predicted by a 6~2 CDM model is far lower than 
the density observed, while the number density for an unbiased CDM model is roughly consistent with 
our value. The COBE measurements of CMBR temperature anisotropies (Smoot et al. 1992) also imply 
that (78 = 1.08 ±0.25 (Efstathiou, Bond & White 1992) and favor an unbiased version of standard CDM. 
The question remains, then, how to reconcile the apparent discrepancy between these observations, and 
observations of small scale structure that require a biased (b ~ 2) version of CDM. 

The currently popular opinion is that standard biased CDM does not quite work (Ostriker 1993; Primack 
1994), and various groups have started to investigate models with an additional free parameter; for example 
CDM models with a non-zero cosmological constant that are spatially flat (Ao/3Hq = 1 — Oo) an d thus 
compatible with inflation (so-called LCDM or A + CDM models; Efstathiou et al. 1990; Kofman et al. 1993; 
Efstathiou 1994), or O = 1 hybrid or mixed dark matter models with varying amounts of both cold and hot 
dark matter (C+HDM ), and with the hot dark matter (HDM) usually assumed to consist of one flavor 
of massive neutrino (Davis, Summers & Schlegel 1992; Taylor & Rowan- Robinson 1992). One popular 
variation of C+HDM that appears to agree with the COBE measurements, and measurements of large scale 
structure and large scale bulk flows assumes an 0=1 universe with ~70% CDM and ~30% HDM. 

But while these popular mixed dark matter models have additional large scale power afforded by an 
HDM component that may help to resolve many of the problems of standard CDM, = 1 versions of these 
hybrid models may be just as vulnerable as standard CDM to the existence of massive clusters of galaxies 
at high redshift. For example, using N-body simulations, Jing & Fang (1994) have calculated the predicted 
abundance of massive clusters at intermediate and high redshifts in such a hybrid model composed of 60% 
CDM (O C dm=0.6), 30% HDM (O H dm=0.3), and 10% baryons (O 6 =0.1), and compared it to the predicted 
abundance for a flat, low-density CDM model with a non-zero cosmological constant (Oo=Ocdm=0.3, 
/i=0.75, and Ao = A/3Hq = 0.7). While both models successfully account for the abundance of clusters 
at low redshift, and agree with the large scale structure, large scale bulk flows, and COBE measurements 
mentioned above, Jing & Fang show that the comoving number density of rich (Coma-like) clusters should 
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Figure 3. Plot of the comoving number density of distant X-ray clusters as a function of X-ray luminosity. 
The three points to the upper left are from Castander et al. (1994) for three z~0.8 optically selected clusters 
observed with ROSAT. The 0.1 — 2.5 keV R OSAT luminosities were converted to 0.3 — 3.5 keV EMSS lumi- 
nosities assuming T~6 keV. The two points to the lower right are for the distant EMSS clusters discussed in 
this paper. The upper point is for all six z > 0.5 EMSS clusters (z = 0.66) and the lower point is for the two 
~z = 0.8 EMSS clusters. For comparison, we also plot the integrated power law luminosity function determined 
by Henry et al. (1992) for their highest redshift shell (0.30 < z <0.60, z = 0.33) of the EMSS clusters. 



evolve rapidly in the hybrid model, decreasing by an order of magnitude from z=0 — ► 0.5 while remaining 
more or less constant for the A + CDM model. By z ~ 0.8, the density of rich clusters in a C+HDM universe 
should be another factor of 10 lower than at z ~ 0.5, and thus a factor of 100 lower that the density at z=0. 
Although the five z^0.8 clusters of Gunn (1990) appear to be in severe conflict with this form of hybrid 
dark matter model, Jing & Fang argue that the density of EMSS clusters at z^0.5 (Henry et al. 1992) 
seems to be consistent with the predictions of C+HDM . One can therefore invoke the usual argument 
that the measured velocity dispersions of the z^0.8 optically-selected clusters may somehow be inflated 
by contamination and projection effects. We note here, however, that the EMSS cluster density as derived 
by Jing & Fang from the data of Henry et al. (1992) did not include the four new z > 0.5 EMSS clusters. 
Moreover, the number density we find at z ~ 0.8 of n z ~ 7xl0 -9 h^ Mpc -3 for the two extremely luminous 
X-ray clusters having L x > 7.5 Xl0 44 erg s _1 may be only marginally consistent with the theory. Discovery 
of additional clusters of this type at these high redshifts will surely present a serious challenge to this variety 
of mixed dark matter model. 

On the other hand, it appears the existence of these high- z clusters does not present a problem for fiat, 
low density A + CDM models, but these models lack a theoretical explanation for a non-zero vacuum energy 
density (see discussion in Kofman et al. 1993). Nevertheless, we also note that a non-zero cosmological 
constant is one of the few ways to adjust the age of a fiat universe so it exceeds the lower bound on globular 
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cluster ages when most modern distance scale methods seem to indicate Ho > 75 (Jacoby et al. 1992; 
Fukugita, Hogan & Peebles 1993; Pierce et al. 1994; Freedman et al. 1994). 
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